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1. Introduction 

Anderson localisation is the term used to describe a generic phenomenon, discovered 
in the late fifties by P. W. Anderson, whereby the addition of a certain amount of disorder 
or randomness in an otherwise deterministic medium causes the waves propagating in 
the medium to become localised in space [2]. For quantum systems, the understanding of 
transport properties requires a thorough study of how the presence of disorder affects the 
nature of quantum states — a question first addressed in an earlier paper of Landauer and 
Helland [10]. Since that time, the one-dimensional case has been discussed extensively 
and has led to a better understanding of the physical mechanisms that are responsible 
for localisation. A remarkable feature of one-dimensional systems is that almost all 
states become localised as soon as there is any disorder. This result, first conjectured 
by Mott and Twose [47] . was made more rigorous by Borland [6] who considered an 
infinite chain of identical localised potentials separated by regions of zero potential. 
Assuming that the lengths of these regions are independent random variables with the 
same probability distribution, Borland studied the growth rate of the wave function on 
a semi-infinite chain with prescribed boundary conditions at one end. He then argued, 
by using a "matching argument", that the positivity of the growth rate implies the 
exponential localisation of the wave functions. A purely mathematical proof using the 
properties of transfer matrices was given by Matsuda and Ishii [44] . 

These physical arguments can in fact be made completely rigorous. The exponential 
growth of the solutions of the Cauchy problem is a crucial feature of the proof that the 
spectrum has no absolutely continuous component [53], and that it is pure-point [31]. 
If ip(x, E) denotes a solution of the Cauchy problem (i.e. a solution of the Schrodinger 
equation on the positive half-line subject to boundary conditions at x = 0) then the 
quantity 

7 (E) := lim (1) 

x->oo X 

is a self-averaging quantity called the Lyapunov exponent of the disordered system. A 
rigorous demonstration of the fact that, under certain hypotheses, 7 also quantifies the 
exponential decay of the eigenf unctions — and therefore that its reciprocal can serve 
as a definition of the localisation length — appears for the first time in Ref. [10J. It 
should be borne in mind, however, that this definition of the localisation length is only 
useful if certain conditions are fulfilled; cases where the definition is inappropriate will 
be considered briefly in £|4] 

In order to illustrate the localisation phenomenon, Anderson made use of a model 
in which the wave function solves a difference equation. The general solution of this 
"tight-binding" model takes the form of a product of random matrices, say, 

IT n := M n M n _a • ■ • M 1 (2) 

where the Mj are independent and identically-distributed square matrices with a 
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common probability measure /i(dM). The quantity 
r E(ln|n n |) 

7 M := Inn 3 

7i— >oo n 

is called the Lyapunov exponent of the product of random matrices and, as we shall 
soon see, is effectively the same as 7. 

Anderson localisation has been a powerful motivation for the study of products 
of random matrices. For this reason, the search for precise conditions on the measure 
ju(dM) that would guarantee the existence and the positivity of the Lyapunov exponent 
7^ has been of particular interest. The main result in this respect is due to Furstenberg 
[28] — a result published the same year as Borland's paper; see also Oseledec [51] whose 
work can be considered as an extension to dynamical systems of the work of Furstenberg 
and Kesten [27] . Further developments of these results, and their application to the 
study of localisation, are described in the works of Ishii [35], Bougerol and Lacroix [8], 
Carmona and Lacroix [11], Lifshits et al. [41 J, Luck [42J, Pastur and Figotin [54J, 
Crisanti et al. [JTJ] and the references therein. 

An important milestone in the development of the theory of localisation was the 
discovery of the relationship between the Lyapunov exponent and the integrated density 
of states N(E). It turns out that the characteristic function 

Q(E) := j(E) - inN(E) , (4) 

viewed as a complex-valued function on R, is analytic in the upper half of the complex 
plane [35] [62] . Interestingly, this analyticity property of the characteristic function was 
in fact exploited much earlier by Dyson in his famous paper on the dynamics of a 
disordered chain [2"T] . 

Products of random matrices appear naturally in several other problems related 
to the physics of disordered systems — see the monographs [TJJJ [42] — and is still the 
subject of active research [24]. The concept of Lyapunov exponent is a very useful 
tool for analysing a large class of systems with quenched disorder, such as magnetic 
systems. A well-known prototype is the Ising model in a random magnetic field, where 
the calculation of the free enegy reduces to analysing an infinite product of 2 x 2 matrices; 
in this context, the Lyapunov exponent is proportional to the free energy per spin. 

The group from which the matrices M n in the product ^ are drawn varies not 
only with the physical context but also with the choice of vector basis. In the analysis 
of the Schrodinger equation in a random potential, a standard choice is to consider the 
vector formed by the wave function and its derivative {ip'{x) , ip{x)). As we shall see 
later on, the evolution of this vector is governed by matrices belonging to the group 
SL(2, K). This is the formulation chosen for example in Refs. [151 ELY] and in the present 
article. If, instead, one chooses to focus on the scattering aspects of the problem then the 
wave function, in a region where V(x) = 0, is a combination of incoming and outgoing 
waves ip(x) = Ae lkx + B e~ lkx . In this setting, it is natural to consider transfer matrices 
T connecting pairs of complex amplitudes (A, B). Current conservation then implies 
that such transfer matrices belong to the group U(l, 1); see the appendix of Ref. [15] 
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for further details, and Ref. [33] for a pedagogical presentation and a review of the 
literature. 

Interestingly, this scattering formulation provides a natural way of generalising the 
one-dimensional model — which is associated with products of 2x2 matrices — to a quasi- 
one-dimensional model where the matrices in the product are 2m x 2m, where m is the 
number of conducting channels. The localisation problem then involves — not just one — 
but rather m (counting multiplicity) Lyapunov exponents 7i ^ • • - ^ 7 m - Since the 
pioneering work of Dorokhov [20J, whose results were later rediscovered independently 
by Mello, Pereyra and Kumar [33], this topic has attracted a lot of attention owing to its 
relevance in the description of weakly disordered metallic wires; see the review [3] . These 
early works relied on the so-called isotropy assumption, namely that each elementary 
slice of disordered metal redistributes the current uniformly amongst the m conducting 
channels. This assumption produces a set of Lyapunov exponents with the behaviour 
7j = 71 [1 + /3 (j — 1)], where (3 G {1, 2, 4} is the Dyson index [3]. The smallest Lyapunov 
exponent, which scales with the number of channels like 71 oc 1/m, is usually interpreted 
as the reciprocal of the localisation length. The ideas of Dorokhov and Mello et al. have 
since been extended to other symmetry classes of disordered models; see the review 
[23]. Another line of research stimulated by these ideas is to look for models where the 
isotropy hypothesis may be relaxed [TJJ 133 13HI EH] , with the aim of studying the passage 
from one-dimensional to higher-dimensional localisation. To close these brief remarks 
on the multichannel case, we point out that the Lyapunov exponent is at the heart of 
numerical studies of localisation that rely on the scaling approach [38] . 

The present paper revisits the interplay between products of random matrices 
of SL(2,M) and one-dimensional Anderson localisation. Whereas most of the 
works cited above emphasise discrete models obtained by making the tight-binding 
approximation [T^l 132] , we shall focus here on one-dimensional continuous models that 
make use of the notion of point scatterer. One familiar example is the Kronig-Penney 
model [39] based on delta-scatterers. We shall see that, by considering a natural 
generalisation of the concept of point scatterer, we arrive at a useful interpretation 
of a general product of random matrices as a model of disorder. We illustrate the 
fruitfulness of this interpretation by exhibiting concrete instances of the probability 
measure fi(dM) of the matrices Mj for which the Lyapunov exponent may be expressed 
in terms of special functions. One of these models is new and has interesting connections 
with Sinai's study of diffusion in a random environment [TJ [59]. Another justification 
for focusing on models involving generalised random point scatterers is that — far from 
being special — they can on the contrary be mapped onto the most general product of 
random matrices in SL(2,R). Hence such models can in principle describe the whole of 
one-dimensional disordered quantum mechanics. 

The remainder of the paper is as follows: ^2] describes the concept of generalised 
point scatterer, explains how models of disorder may be constructed from them, and how 
completely general products of matrices may be associated with such models. We then 
list some known solvable cases. £j3] is devoted to the analysis of a new example where 
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the Lyapunov exponent may be expressed in terms of the generalised hypergeometric 
function. The paper ends in £j4] with a discussion of the limitations of the Lyapunov 
exponent as a measure of localisation in disordered systems. 



2. Generalised point scatterers and products of random matrices 

2.1. Point scatterers and transfer matrices 

A point scatterer is the idealised limit of a potential whose action is highly localised. 
The most familiar example is the delta-scatterer at a point, say Xj. In the context of 
the Schrodinger equation 

-i)"{x) + V(x)xp(x) = Exp(x), x>0, (5) 

this potential vanishes for every x ^ Xj and is defined at x — Xj via the boundary 
condition 



ip'ix^ 

1p(Xj+) 




(6) 



where Uj is the strength of the scatterer, i.e. V(x) = ^2jUj5(x — Xj). This delta- 
scatterer is generalised by replacing the 2x2 matrix on the right-hand side by an 
arbitrary matrix, say Bj. Thus 

Conservation of the probability current requires Bj e SL(2,R). If, instead of a single 
scatterer, we consider a sequence of scatterers placed at the points 

X\ < X2 < ■ ■ ■ 

then, assuming E = k 2 > 0, the solution of the Schrodinger equation may be expressed 
in the form 



ip'(x n +) 

where IT n is the product ([2j with 




(8) 



sin 9j cos 9j 




and 



$ j ft' ^3 ' ^3 ' "^'3 I 1 *^ 3 ' 

Now, by applying the Gram-Schmidt algorithm to the columns, every 2x2 matrix M 
with unit determinant may be expressed in the form 

, r , cosfl -sin^ \ / e w , , . . ,. , 

M= l sin e cose II n - II n , I < 10 ) 
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for some real parameters 9, w and u. This is the Iwasawa decomposition of SL(2, M.) into 
compact, Abelian and nilpotent subgroups. Therefore, by working with point scatterers 
of the form 

o W i «, \ 

3 I o <r w i j I o l J y ' 

and taking E = k 2 = 1, we obtain a correspondence between general products of 
matrices IT ra in SL(2,R) and the Schrodinger equation. The parameter k may be easily 
reintroduced by simple dimensional analysis. 



2.2. Product of random matrices — The Riccati variable 

There are several options in defining a disordered quantum system with point 
scatterers. The randomness may be in the strength of the scatterer, i.e. in the matrix 
Bj, or in the position of the scatterer, i.e. in the coordinate Xj or in both the strength 
and the position; see [JT] where various models are reviewed. 

We shall henceforth confine our attention to the particular case where the spacing 
between consecutive scatterers, i.e. the angle 9j = k£j, is exponentially distributed: 



F(£j eS)=p / e~ p "d£ (12) 

where 1/p > is the mean spacing. This is the situation that arises when impurities 
are dropped uniformly on R with a mean density p. 
The equation satisfied by the Riccati variable 



IS 



and 



Z'(x) = —E — Z 2 (x) for x i { Xj } (14) 



Z(xj+) = Bj {Z(xj-)) for j EN (15) 



where Bj is the linear fractional transformation associated with the matrix Bj, i.e. 

Bj{Z) :=e 2w i {Z + Uj) . (16) 

Because the spacing between consecutive scatterers is exponentially distributed, Z is 
a Markov process, and it was shown by Frisch and Lloyd [25] (see also [15]) that the 
stationary density f(Z) satisfies 

d 



dZ 



[(Z 2 + E)f(Z)} +p I ii B (dB) 

JSU2,R) 



f{B-\Z))^^-f[Z) 



(17) 



'SL(2,1 

where fiB is the probability measure of the random matrix B. 

The relationship between the Riccati variable and the Lyapunov exponent is easy 
to establish — at least heuristically. A completely rigorous treatment would follow the 
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lines of Kotani's work [37J. Given the particular form of the upper-triangular matrix 
Bj, we have 

In \ip(x n +)\ = -w n + In \tp(x n -)\ 

= ~w n + / Z(x) dx + In \ip(x n -i+) | 



71 fXn 

= --- = -^2wj+ Z(x)dx + \n\ip(x 1 -)\ . (18) 
Dividing by x n and letting n — » oo, we obtain 

lim ln l^( x ")l = _ lim 1_ y + lim ! \ Xn Z m dx . (19 ) 



Since 



n->oo £ n->-oo £ A — ' n->oo X 

n-l 



where the £j are independent and identically distributed with mean 1/p, the law of large 
numbers implies that 

x n ~ n/p almost surely, as n — > oo . (20) 

Then, by the ergodic theorem, 

j=-pE(w)+ [ Zf{Z)dZ (21) 

almost surely, where / is the density of the stationary distribution of the Riccati variable. 
The upshot is that the Lyapunov exponent of the system may in principle be computed 



by solving the Frisch-Lloyd equation ( 17) for the stationary density, and then performing 



the integral in Formula (21). 



2.3. Correspondence with the Furstenberg theory 

It is instructive to compare this method of calculation with that based on 
Furstenberg's theory [261 EQ]- The concept of direction, through the projective space 
P (R 2 ), plays a prominent part in that theory. This is due to the fact that the product of 
matrices grows if and only if the columns of the product tend to align along a common 
direction. In M 2 , a direction may be parametrised by a single number; let us choose the 
reciprocal zgKU {oo} of the slope. The columns of the product n n have directions 
that are random. When a vector of random direction, say z, is multiplied by a random 
matrix 

a new vector is produced, whose direction is the random variable 

MM == ^ • (23) 
cz + a 
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We say that the distribution, say z^(dz), of z is invariant under the action of matrices 
drawn from \x if the distributions of the old and the new directions are identical, i.e. 

M(z) l =z. (24) 

In particular, if i^(dz) has a density, say f/j,(z), then this equality in law translates into 
the following equation for / M : 

dA^ 1 , s 



SL(2,1 



MdM) {f.oM- 1 ) (z) 



dz 



(25) 



where, with some abuse of notation, z is no longer random. Knowing f^, the Lyapunov 
exponent 7^ of the product U n may then be computed by the formula 



7m 



fdzfrWf //(dM)ln 

JR >/SL(2,K) 





(26) 



The calculation of 7^ when M is of the form ( 10 ) and 9 is exponentially distributed 
with mean 1/p, independent of w and u, may then be related to the calculation of 7 in 
the previous subsection as follows: 

The projective variable z is the Riccati variable Z with E = k 2 = 1. 



m 



The Furstenberg equation (25) for the invariant density reduces to the Frisch-Lloyd 



equation (17) for the stationary density. 



The Furstenberg formula (26) for 7 M reduces to Formula (21) for 7 and 
1 



7m 



P 



■7- 



(27) 



The factor 1/p is due to the difference between the definitions: whereas ([T]) involves 
the distance variable x, the growth of the matrix product is measured with respect 
to the index n in pi) — x and n being related by ( 20 ) . 



2.4- Known solvable cases involving random point- scatter ers 



There is no systematic method for solving the integro-differential equation (17), but 



some solvable cases have been found and are listed in Table Q] All but the first of these 
cases have in common that one or both the parameters w and u in the expression for 
the matrix B are exponentially or gamma distributed. The density of the exponential 
and gamma distributions satisfies a differential equation with constant coefficients, and 



this leads to a trick for reducing (17) to a purely differential form. This trick will be 
illustrated in the next section. 

Here, we consider the simplest example and merely write down the result. This 
example consists of taking 



B 




(28) 
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where u is exponentially distributed with mean 1/q. Then the stationary density solves 
A [(^ + k 2 ) f(z)] - pf(z) + q(z 2 + k 2 ) f(z) = qN . (29) 

The constant A" = N(E) appearing on the right-hand side is the integrated density of 
states per unit length of the quantum model; the Rice formula 

N(E) = lim z 2 f(z) (30) 

\z\— >OD 

expresses its relationship to the tail of the stationary density. After integration, we 
obtain 



/(*) 



qN 



cxp 



p z 
—qz + — arctan — 
k k 



exp 



p t 
qt — — arctan — 
k k 



dt. 



z 2 + k 2 

The integrated density of states is then determined by the requirement that / be a 
probability density. The final formulae for A" and 7 may be expressed neatly via the 
characteristic function Q of the disordered system: for E = k 2 > 0, 



Q(E) = 2ik 



ip 1 
2k '2 



-2\kq) 



W 



-ip 1 (— 2ikq) 

2k '2 



(31) 



where W a ,p is a Whittaker function [32]. This formula, which was discovered 
by Nieuwenhuizen [50] using a different approach, corresponds here to the case of 
"repulsive" scatterers (u > 0); it is easily adapted to the "attractive" case u < 
by simply changing the sign of q in Eq. (31) [15J. 



Density of I Density of u Density of w Reference Special function 



qe -qu | 
-,2 



-1/p) 

pe-P e l M+ {£) 

pe- pe l R+ (i) q 2 ue-i'l 

pe-P e l R+ (£) 8(u) 



M.AU) 



5{w) 
5{w) 
.(«) 6(w) 

qe' qw l 

q e -2q\w\ 



[501 H5] Whittaker 

[501 US] Whittaker 

[15] Hypergeometric 

this article generalized Hypergeometric 



Table 1. Known solvable cases involving random point scatterers. The notation 1^4 (x) 
is used for the function that equals 1 if x £ A and otherwise. 



Besides the Lyapunov exponent, the low-energy behaviour of the density of states 
is also of interest; it has been analysed for point scatterers with random positions in the 
following cases: {%) when w = and the probability density of u is supported on the 
positive half-line but is otherwise arbitrary [371 HI] : (ii) in the converse situation where 
u = and it is w > that is random [16J. 

2.5. Continuum limit and models involving Gaussian white noises or Levy noises 

No exact solution is known in the case where all three parameters i = 8, w and u 
are random. Some insight into the general case may however be gained by considering 
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the continuum limit of an infinitely dense set of point scatterers with vanishingly weak 
strengths, i.e. 

6 , w , u — > . 

This limit was studied very recently in Ref. p2] in the most general case where the three 
parameters of the Iwasawa decomposition ( 10 ) are correlated — the matrices Mj being 
still mutually independent. This recent work corresponds to a generalisation of the 
study of the Schrodinger equation with a potential V(x) = r](x) + W(x) 2 — W'(x) that 
combines two (possibly correlated) Gaussian white noises rj and W [33]; see Ref. [T7] for 
a detailed discussion. 

Finally let us mention a generalisation of the model in another direction: by letting 
the density of scatterers tend to infinity, it is possible to study certain singular limiting 
measures for the strength of the scatterers. Such cases may be described in terms of 
Levy noises and were studied in 



3. Supersymmetry: a new solvable case 



In this section, we shall restrict our attention to the particular case where the point 
scatterers are described by transfer matrices of the form 











As shown in Ref. [To] , this corresponds to considering the potential 

V = w 2 - w 

where W is a superposition of delta-scatterers 



(32) 



(33) 



A potential of the form (33) is said to be supersymmetric with superpotential W . Hence 
we shall refer to the point scatterer associated with (32) as "supersymmetric". This case 



has many interesting features: the corresponding Hamiltonian is factorisable and, when 
the superpotential is deterministic and has the so-called "shape-invariance" property, 
the discrete spectrum may be obtained exactly [IS] . 

The disordered case is of interest in several physical contexts [13] : in particular, the 
Schrodinger equation can be mapped onto the Dirac equation with a random mass — 
a model relevant in several contexts of condensed matter physics; see the reviews 
[HI [2S1 SI] and the introduction in Ref. [HI]. The model is also closely connected 
to diffusion in a random environment [7]; indeed, if we set 



ip(x) = exp 



W(x) dx 



U{x) 



then 



U"{x) - W(x) U'(x) = X U(x) with A 



E 
2~ 



(34) 
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When E < 0, this is the equation satisfied by the Laplace transform of a hitting time of 
the diffusion in the environment J W(x) dx; see for instance [9]. 



After integrating with respect to z, the Frisch-Lloyd equation ( 17) for the stationary 
density becomes 



(z 2 + E)f(z)+p / dwg(w) / f(t)dt 



N(E) (35) 



where g is the density of the random variable w. Our strategy for computing the 
Lyapunov exponent will make use of three simple observations: firstly, the equivalence 



between the supersymmetric Schrodinger equation and Equation (34) implies that the 
spectrum is necessarily contained in IR + . It follows that N(E) = for E < and 
that the stationary density / is supported on the positive half-line; this facilitates the 
calculation of the Lyapunov exponent. The second observation is that the characteristic 
function (|4]) is an analytic function of E, except for a branch cut along the positive 
half-axis. Therefore, if we find an analytical formula for the Lyapunov exponent 
when E < 0, analytic continuation will furnish a formula for the case E > 0. The 
third observation — which we alluded to earlier and whose exploitation goes back to 
the works of Nieuwenhuizen [50] and Gjessing and Paulsen [28j — is that the integral 



term appearing in (35) may be eliminated if g solves a linear differential equation with 



(piecewise) constant coefficients. 

With these points in mind, let us set 



and take 



E = —k < 



g(w) =qe- 2qlwl . (36) 



If we write f p , q ,k to indicate the dependence of the stationary density / on the parameters 
p, q and k, then we have the following identity: 

• (37) 

It is therefore sufficient to consider the case k — 1, i.e. E = — 1. It is then easily deduced 
from the evenness of g that 

f(z) = -J(l/z). (38) 

z A 

As a consequence, to know f(z) for z ^ 1 is to know f(z) for < z < 1, and vice-versa; 
we shall at times implicitly make use of this helpful property. 

3.1. Reduction to a differential equation 

Following the recipe outlined in [15], we shall now deduce from the Frisch-Lloyd 



equation (35) a differential equation for the density /. To this end, we introduce the 
kernel 

K(y) = -sign(y) V 2 ^l (39) 
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and note, for future reference, that K satisfies 

K'(y) = - S ign(y)2qK(y), y^O, 

subject to 

K{0±) = t\. 

After permuting the order of integration in the Frisch-Lloyd equation, we find 



Set 



z 2 -l)f(z)+p / K -In- /(t)dt = 



<p{z) := (z 2 - l)f(z) 



12 
(40) 
(41) 

(42) 



and write 



f K (i ln I ) /(() * = I A ' S 9 /(t) * + f K G ln S /(() * • 



Differentiation of (42) with respect to 2 then yields 

9 



<p'(z)+pK(0+)f(z)-p- 



z 







A I 1 h/ ) /(/ )d/ 



-pif(0-)/( Z )+p 



9 



A I iln|] /f-/-)d/- = 0. 



(43) 



Hence, in view of (41), 



[^(*) - p f(z)] - pq Z K Q ln ^ /(t) dt + pq J™ K^-\n |) /(t) dt = 



Another differentiation leads, after use of (41), to 
d 



-{z y{z)-pf{z)]}+p q - 



K 



\hx- t ) f(t)dt . 



(44) 



The integral term may then be eliminated by making use of Equation (42); the result 
is the second-order linear differential equation 

d 



dz 



z \ (p -p 



-<^V = o. 



(45) 



3.2. The case q — 1 

It may be verified by direct substitution that, when q = 1, the differential equation 



(45), expressed in terms of the unknown /, has the solution 

1 + pz + z 2 



f P ,iA z ) 



A 1 

p — 2 z 



1 - 



z-lMOz + l) 1 ^ 



(46) 
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1 - z \ 5 l+pz + z 2 
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2A 
J () 

2A 



1 



l + z 



1-z 2 



p-2 



1 — x 2 
1 — x 2 



P\ £-1 

- X 2 

2 



dz 

z 

dx 



2A f Z" 1 1 - ar§ -1 2 / p 

= ^i 2 y -T^ dx+ p-v-2 

The value of this last integral is given explicitly by 
function Hence 

1 l#(p/4) -#(1/2) 2 



(47) 



(4* 



3.269, in terms of the digamma 

(49) 



A 2 p/4-1/2 
The case p = 2 may be studied by letting p — > 2 in Equations ( 46p9 ). We get 



f2,i,i{z)=A 



2z 



In 



z + l 



z-1 



z + l) 



and 



1 



71 



A 4 



- 1 . 



(50) 



Plots of the density (46) for several values of p are shown in Figure [TJ The integral on 




Figure 1. Plots of the stationary probability density f(z) for q = 1 and energy E = — 1. 



the right-hand side of Formula (21) may also be computed exactly; the end result is 

.2x P 



1 



2A- 



(51) 



2 *(i)-*(!) ■ „ 

Analytic continuation to positive energies E = k 2 > consists of replacing k by — ik; 
this yields 



p 



2 *(l)-*Q)-?-i' 

Various limits may be analysed with the help of these expressions. Using 



(52) 



z^o z 6 
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0.577. 

n(-k 2 ) 



0{k 



(53) 



where C = 0.577. . . is the Euler-Mascheroni constant, we deduce from (|51|) 

^ i. 

Analytic continuation then shows that N(k 2 ) ~ k/ir in the limit k — > oo, as expected 
from the free case. Also, the high-energy Lyapunov exponent tends to a constant: 
7(fc 2 ) ~ p (In 2 - 1/2); see Fig. [2j 

In the low-energy limit, by using 

$(z) = In z — ^ — h 

we obtain 

p/2 



0{z 



n(-k 2 ) 



fc -,o in(f ; 



c 



+ 



In 



(54) 



Analytic continuation to positive energies allows one to recover the Dyson singularity 
of the integrated density of states and the corresponding vanishing of the Lyapunov 
exponent: 



N(k 2 



9 



and 7(/c z 



9 



2\n 2 k ,v " ' |lnfc| 

with g = p/2. In this k — >■ limit, we thus recover as expected the same results as 
in the case where the superpotential W is a Gaussian white noise; the characteristic 
function for this case is recalled below. The fact that 7 — > 0+ as E — > 0+ means that 
the eigenstates become delocalised (or extended) in that limit; see the discussion in 
Section |4} Plots of 7 and iV for positive energy are shown in Figure [2} 



rg 0.20 




point scatterers 
for p =2 



g.w.n. 




p=2 



Figure 2. Plots of N and 7 against k = yE for p = 2 and q — 1. The 

integrated density of states exhibits the Dyson singularity. For comparison, the results 
corresponding to the case where W is a Gaussian white noise of strength g = p/2q 2 = 1 
are also shown as blue dashed curves. The thin line is the asymptote of the Lyapunov 
exponent. 
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3. 3. The characteristic function for arbitrary q > 

Although the calculations are somewhat tedious, the more general case q > is 
also tractable. The results are summarised in the following 

Theorem 1. The invariant probability density is given by 

z Jo v x 
where C is a normalisation constant and M a ^ is a Whittaker function. 

Corollary 3.1. For E = k 2 > 0. 



n (B + ,o +) = -tti±i ^(^-S^ + ^ + ^+l- 1 ) 

9 3 *i(g + £- §,9 + 1,9 + 1; 2« + l, g + §;l) 



Proof. See |Appendix A 



□ 



Although this formula for Q in terms of a known special function is pleasing, the task 
of extracting from it concrete information is not entirely straightforward. We proceed 
to discuss various limits. 

The integrated density of states behaves in the limit E — > +oo as in the free case. 
We may obtain the behaviour of 7 by using the perturbative approach described in 
Ref . [1] . When the superpotential W is given by a superposition of delta-functions, the 
high energy Lyapunov exponent of the supersymmetric Hamiltonian tends to a finite 
limit given by 700 = pE(lncoshw). Therefore we obtain 



l{E) 



P 
2 



q 



-*(? + i 



q + l 



as E — > 00 . 



(55) 



For example, in the case q — 1, we recover the limit 7 ~ p (In 2 — 1/2) discussed earlier. 

The low-energy limit may be conveniently studied by using the integral 
representation of the characteristic function 

dx 

Fx 



tt(-k 2 ) = k 



dx 



(56) 



M^ q (x)K (x/2) 



following from (A. 25). We set E 



-1 for simplicity. The p — > 00 limit provides 



the E — > behaviour, after reintroducing k. We now analyse the integrals in this limit. 
There holds [32] 

(57) 



M| i9 (x) ~ x q+1/2 as x ->■ 

and, since the Whittaker function solves the differential equation 

x 2 —px + (4q 2 — 1) 



w"(x) 



Ax 2 



w(x) 



we expect the asymptotic form (57) to hold for x smaller than (4g 2 — l)/p (this reasoning 
assumes that 4q 2 — 1 > 0). Moreover, we see from the differential equation that 
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the Whittaker function is a rapidly oscillating function in the interval where 
x± — | ± a/(| ) 2 — 4q 2 + 1; in the limit of interest here we have x_ ~ (4g 2 — l)/p and 
x + ~ p. Thus we expect that the contribution to the integrals of the interval [ X— , X-f-] is 
negligible. For x exceeding x + we have 
1 



x 



MeJx)K s (x/2) 



-2-1 
~ X 4 



from which we get 
(58) 



Finally both integrals in (56) are dominated by the interval [0,x 

n(-i) ~ i±l 2 

p->oo g x_ln(l/x_) hip 

Re-introducing k, we indeed get the expected behaviour k 2 ) ~ — 1/ln/c for fc — > 0. 

The limit p — > oo and g — > oo keeping p/q 2 fixed corresponds to a case where W(x) 
converges to a Gaussian white noise of strength g = pE,(w 2 ) = p/2q 2 — a model that has 
been solved in Ref. [7J E2] (see also |4"T]). Therefore we expect that the characteristic 



(59) 



function (56), i.e. Eq. (A. 30), converges in this limit towards the expression 

Kxik/g) 



n(-k 2 ) 



k 



p— >oo , q= 



/p/2g 



K (k/g) 



4. Concluding remarks: limitations of the Lyapunov exponent for the study 
of localisation 

To close this short review, we point out the limitations of the Lyapunov exponent 
as a means of characterising the localised nature of the spectrum in a disordered system. 
A first observation is that the definition ([T]) relies on the assumption of a linear increase 
of In \ip{x] E)\, i.e. on the application of the central limit theorem. It is however 
possible to consider models of disorder such that In \if)(x\ E) \ scales as x a , leading to the 
phenomenon of sub-localisation (for a < 1) or super- localisation (for a > 1); this has 
been studied in, for example, Refs. [U El H31 [58] (see also references therein). However, 
this limitation is easily overcome by extending the notion of the Lyapunov exponent so 
that it characterises the typical length over which eigenstates are localised. 

Another, more serious, limitation of the usefulness of the Lyapunov exponent 
pertains to the possible existence of extended states in one dimension. This might occur 
when the disordered potential is correlated [36] — as happens for instance in the dimer 
model. It might occur also for symmetry reasons, as is the case for the supersymmetric 
Hamiltonians studied in Refs. H3J HU [331 El] and in section [3] there, the Lyapunov 
exponent j(E) decays to zero in the limit of zero energy. In this respect, it must 
first be pointed out that the vanishing of the Lyapunov exponent is a necessary but 
not a sufficient condition for delocalisation; sometimes it only signals sub-localisation. 
Secondly, several works on disordered supersymmetric quantum mechanics have shown 
that the Lyapunov exponent, when it vanishes at the bottom of the spectrum does not 
provide the relevant length scale characterising the localisation of eigenstates at low 
non-zero energies. As mentioned already in section [3j when E(W) = 0, the Lyapunov 
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exponent vanishes like j(E) ~ g/\ lnE\ as E — > [13J, [H]. This suggests that low- 
energy eigenstates are localised on a scale commensurate with | In E\/g. However, the 
study of the ordered statistics problem [60l El] , as well as the analysis of other physical 
quantities such as the averaged Green's function at non coinciding points [291 E] or the 
boundary-sensitive averaged density of states [51], show that eigenstates are localised 
on a much larger scale, namely | In E\ 2 /g. 

The inability of the Lyapunov exponent to capture the localisation property for 
energies E ~ E c where ^y(E c ) = may be related to the fact that it is defined, via 
Eq. 0, in terms of a solution ip(x; E) of the Cauchy problem that vanishes at just one 
boundary. This definition may indeed fail to describe the interesting properties of the 
real eigenstates — which are solutions of the Schrodinger equation vanishing at the two 
boundaries x = and these eigenstates become less localised. 
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Appendix A. Proofs of Theorem [T] and its corollary 



Let us look for a solution of the differential equation ( 45 ) of the form 

p 



Then 



y" + 



l + z 



2g + l + _p/2_ + -p/2 
z z — 1 z + l 



y' + 



pq 



z(z-l)(z + l) 



y = 0. 



(a.i; 



The solution we require has various representations in terms of special functions. 



Appendix A.I. Heun's function 

Heun's differential equation is [5 
~V 



y + 



5 

+ 7 

Z Z — 1 



+ 



z — a 



y' + 



a (3z — q 



z(z — l)(z — a) 



y = 



(A.2) 



where the parameters satisfy the relation 

a + + l = n + 5 + e. 
Heun's function is the particular solution analytic inside the unit disk: 



Hl(a,q;a,/3,r),5;z) := 



y n z 



\z\ < 1 



(A.3) 



n=0 
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where the y n satisfy the recurrence relation 
a(n + 2)(n + 1 + r])Vn+2 

- [(n + \){n + ?7 + 5) a + (n + l)(n + i] + e) + q] y n+1 

+ (n + a)(n + fi)y n = (A.4) 

with the starting values 

y = l and ^ = — . (A. 5) 

ai] 



Comparing Equation ( |A.l ) with Heun's equation (A. 2), we find the particular solution 



y = Hl(-l,-pq;0,2q,2q + l,p/2;z). 
The recurrence relation for the coefficients y n of this solution is 

(n + 2)(n + 2q + 2) y n+2 = p(n + q + 1) y n+1 + n(n + 2q) y n (A.6) 
with the starting values 

y = 1 and y 1 = - . (A.7) 

It follows in particular that y is strictly positive for z > 0. Hence, for z G [0, 1), 

f ^ z) = Y^\T+~z) Hl (-h-pq;0,2q,2q+l,p/2;z) (A.8) 



where C is a normalisation constant. The property (38) provides the obvious formula 
in the interval z > 1. 

Appendix A. 2. Gauss' hypergeometric function 
Set 

a := | + q - - , b := | and c := g + - . 

Then, by comparing the Taylor series, we find that 

Hl(-l,-pq;0,2q,2q + l,p/2;z) 
(1 + z) 2b 

= —, ! —{{a + l)(b+l)z 2 {l-z 2 ) 2 F 1 (a + 2,b + 2;c + 2;z 2 ) 

c{c +1) 

+(c + 1) [c(l - 2) - 6z] (1 + z) 2F1 (a + 1, b + 1; c + 1; z 2 ) } . (A.9) 
The contiguity relation 

t(l -t)(a + 1)(6 + 1) 2 Fx(a + 2, 6 + 2; c + 2; t) 

+ [c - (a + 6 + l)t] (c + 1) 2 Fx(a + 1, 6 + 1; c + 1; t) 

-c(c + 1)2^1 (a, 6; c; t) = (A.10) 
leads to the simpler expression 

Hl(-1, -pq; 0, 2q, 2q + l,p/2; z) (A.ll) 

= (1 + z) 26 | 2 Fi (a, 6; c; z 2 ) - h -z{\ - z) 2 F X (a + 1, 6 + 1; c + 1; z 2 ) J . 
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then yields 



2 F[(a, b] c; t) = — 2 F 1 (a + 1, b + 1; c + 1; t) 

c 



Hl(-l,-pq;0,2q,2q+l,p/2;z) 
= (1 + z) 2b \ 2 F 1 (a, b; c; z 2 ) - -(1 - z) 2 F[ (a, b; c; ; 

= (1 + z) 2b | 2 F! (a, 6; c; z 2 ) - ^(1 - 2) (a, b; c; z 2 ) 

2a(l - zf a - x 2 F 1 (a, 6; c; z 2 ) 
1 2 Fi (a, 6; c; z 2 ) | . 



-1 (1 + ^ 
2a (1 - z; 



,26 



2a-l 



\2a 



+ (!-«) ^ 



So we obtain the "compact" representation 

Hl(-1, -pq; 0, 2g, 2g+l,p/2; z) = -(l-z) 2 ^ 

where 



1 + z 
l-z 
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(A.12) 
(A.13) 



(l-z) 



§ + 2g-l 



V 



1 p 



1 



2'4' 



(A. 14) 



i?»,(A.15) 



(A.16) 



Appendix A. 3. The associated Legendre function 

Formula (6) in [22], §3.13, expresses the hypergeometric function in terms of the 
associated Legendre function: 

r (q + Q|_7 (coshr?) = 2 v ^e i (?^)- r f | + q + ^ 



x 



-v 1 9 



(1-e-^ 



1 + e' 
1 - e-^ 



. p 1 

V\ 4~2 



^ (e^) 



(A.17) 



Equation (8) in |32j, §7.621, gives a useful integral representation of the Legendre 
function in terms of a Whittaker function: 

dx 



-U*+V*)M P _x a _i(x) 



4 2 '1 2 



r(2g) 



2 _, ^ e^S-l' 



T(g + 1-|: 



1 + zJ ^ 



Then, by Equation (A.17), we obtain 



(l-z) 



q l-z 

F'(z) 



z + l/z 



r(g+l/2) 



(A.18) 



1 + z 



4A(I-| + ?) r(2g) 



x 



e * 



(z+1/2) M E _i„ _i(x) 



4 2 ,y 2 



xz 



1 - z 

2z 



dx . 



(A.19) 
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Now, 



-4(s) 



1 - z 
2z 



Mp_i „ i(x) 



A 9 '*? o 



+ g-a„ , , 1 



1 + z< 



q 1\ 1 1 + 

x 2/2 4^ 2 



M E _i„ _i(x) (A.20) 



Iq^fx z 4 z 4 2 2 4,2 2 

where we have used Formulae (13.4.28) and (13.4.32) from [1J to obtain the last equality. 
The proof of Theorem [l] follows after reporting this in Equation (A. 19) and using 
integration by parts for the term involving the derivative of the Whittaker function. 



Remark. Formula (5) in [22], §2.8, expresses F p ^ in simple terms: 

-l 

+ 1 



1 



V 



1 - z 



1 + z 



(A.21) 



It then follows from Equation (A. 17) and from the recurrence relation for the associated 



Legendre function (see Formula (18) in [22], §3.8) that, for q e M, F p>q may be expressed 
in terms of elementary functions. For example, we have 

1/8 (1-z) 2 



and 



F P Az) 



(!-§)(!+§) - 

3/32 



l + z 



(4 2) (4 2) (4 2) (4 2) 



'l-z) 



x 



l + z 



1 2 + 1 



Formulae for / p ,i,i and / p ,2,i are easily deduced. 



1 U + i 



(A.22) 



. (A.23) 



Appendix A. 4- The characteristic function 
Define 

POO 

fp,q,k( S ) '■= Z~ S f p ^ k (z)dz. 

Jo 

Using Theorem [I] and changing the order of integration, we find 

''''' 2 AC f°° MUx)W , s (x) d ' 



fp,q,l( S ) — 2 C 



X 



X 



Next, Formula (7) in [22], §6.9, says 



Hence 



/ Mil (s)=2^C(-l) 



M_ P q (-x)W 0>s (x) 



dx 



(A.24) 



(A.25) 



(A.26) 



The integral on the right-hand side is of the same form as that in [32], §7.625, Formula 
(1) — albeit outside the range indicated since in our case a = —1 and j3 — 1. 



ID Anderson localisation and products of random matrices 



21 



Nevertheless, direct numerical verification indicates that the formula does remain valid 
and so 

T(q + l + s)T(q + l- s) 



X3F2 I q 



C 



r(g + §) 



1 p 3 

^ - ^, Q + ^- + s, q + l - s; 2q + 1, q + 1 



The normalisation condition 

LM = 1 

yields 

C 



r(? + §) 



2v^r(g + l)2 3 F 2 ( g +i-f,g+l,g + l; 2g+l,g+§; l) 

and so we deduce the formula: 

f 1 \ _ T(q + l + s)T(g + l-s) 
h,iA s ) r(g + l) 2 

3*2 (q + I - I , q + 1 + s, q + 1 - s; 2g + 1, g + §; l) 



x 



$J F 2 (9+§-f,?+l,g+l; 2g+l,g+§; l) 



(A.27) 



(A.28) 



(A.29) 



for — 1 < Res < 1. Then, using Formula (21 ) for the Lyapunov exponent, together with 

(A.30) 



the identity (37) and the fact that N(E) = for E < 0, we eventually find 



Vt(-k 2 ) = k 



q + 1 3F 2 (q+l-^,q + 2,q;2q+l,q+l;l) 



^,g+l,g+l; 2g+l,g+|; l) ' 



1 3*2 (q+ \ 4fe , 

The analytic continuation from negative to positive energy consists of replacing k by 



—ik. This completes the proof of Corollary 3.1 
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